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ABSTRACT 

We have used the InfraRed Array Camera (IRAC) onboard the Spitzer Space 
telescope to search for low mass companions of the nearby debris disk star e Eri- 
dani. The star was observed in two epochs 39 days apart, with different focal 
plane rotation to allow the subtraction of the instrumental Point Spread Func- 
tion, achieving a maximum sensitivity of 0.01 MJy/sr at 3.6 and 4.5 /im, and 
0.05 MJy/sr at 5.8 and 8.0 /j,m. This sensitivity is not sufficient to directly de- 
tect scattered or thermal radiation from the e Eridani debris disk. It is however 
sufficient to allow the detection of Jovian planets with mass as low as 1 Mj in 
the IRAC 4.5 /xm band. In this band, we detected over 460 sources within the 
5.70' field of view of our images. To test if any of these sources could be a low 
mass companion to e Eridani, we have compared their colors and magnitudes 
with models and photometry of low mass objects. Of the sources detected in at 
least two IRAC bands, none fall into the range of mid-IR color and luminosity 
expected for cool, 1 Gyr substellar and planetary mass companions of e Eridani, 
as determined by both models and observations of field M, L and T dwarf. We 
identify three new sources which have detections at 4.5 [im only, the lower limit 
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placed on their [3.6]-[4.5] color consistent with models of planetary mass objects. 
Their nature cannot be established with the currently available data and a new 
observation at a later epoch will be needed to measure their proper motion, in 
order to determine if they are physically associated to e Eridani. 

Subject headings: circumstellar matter — infrared: stars — planetary system - 
stars: individual (e Eridani) 

1. Introduction 

e Eridani (HD 22049, GJ 144, IRAS 03305-0937, HIP 16537), along with Vega, Fo- 
malhaut and (5 Pictoris, is one of the "Fabulous Four" nearby main sequence stars found 
by IRAS to have a cool infrared excess (Gillett 1986). This excess is explained as thermal 
emission from cool dusty material. A careful analysis of the IRAS data has resolved the 
emitting regions for some of these stars, showing a disk-like structure with a typical size of 
a few hundreds AU, analogous to the Kuiper belt of our own solar system (Beichman 1987), 
but with a total mass of dust several orders of magnitude larger. Given that the timescales 
for the removal of primordial material (leftover from the star formation process) is much 
shorter than the age of these stars, it is believed that these structures - debris disks - are 
generated by the collisions between planetesimals and by cometary activity. As such, debris 
disks are an indirect evidence of planetary system formation (Backman & Paresce 1993). 

With an age below 1 Gyr (850 Myr according to the most recent estimate by Di Folco et 
al. 2004, based on fitting VLTI angular diameter measurements with evolutionary models) 
and a K2 V spectral type, e Eridani is a relatively young star with a mass slightly lower than 
the Sun. The proximity of the system, 3.22 pc distant from the Sun (Perryman et al. 1997), 
offers the rare chance to study in detail the final phases of planetary system formation for a 
star similar to the Sun. 

The debris disk of e Eridani was first imaged at 850 ^m by Greaves et al. (1998) using 
the Sub-millimeter Common-User Bolometer Array (SCUBA) at the James Clerk Maxwell 
Telescope in Hawaii. The SCUBA images show a nearly face-on ring structure, extending as 
far as 35" from the star, with a maximum of emission located at a radius of ~ 18". At the 
e Eridani distance this corresponds to a disk radius of over 100 AU, peaked at the distance 
of 60 AU from the star. The nearly face-on orientation of the disk is in agreement with 
measured projected rotational velocities (v sin i) indicating that the star is seen nearly pole- 
on (Saar & Osten 1997), with an angle of i ~ 30° ± 15° for the pole of the star (deprojected 
adopting the rotational period of 11.68 d from Donahue et al. 1996). 
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The 1998 sub-millimeter map of e Eridani has revealed the presence of a few clumps of 
emission along the ring. Subsequent maps obtained in 2002 with SCUBA at 850 and 450 /im 
(Greaves et al. 2005), and at 350 //m with the SHARC II camera at CSO (Wilner & Dowell, 
personal communication), have confirmed the presence of some of these clumpy structures. 
Greaves et al. (2005) attempted to measure the orbital motion of the structures detected 
in their 1998 and 2002 SCUBA images, and preliminary results suggest a rotation of ~ 6° 
counterclockwise (Greaves et al. 2005). The clumps have been interpreted as evidences for 
the presence of a sub-stellar companion in resonant interaction with the dust: if confirmed, 
their measured motion would require an orbital radius of the resonant body of ~ 20 AU. 
Sophisticated dynamic models of the sub-millimeter ring suggest an orbital radius for a 
Jupiter-mass planet between ~ 40 AU (Ozernoy et al. 2000) and ~ 60 AU (Quillen & 
Thorndike 2002) from e Eridani (12 - 20" at e Eridani's distance) to have the required 
perturbing effect. 

Radial velocity measurements have indeed discovered a Jovian class planet (Msini = 
0.86 Mj 1 ) with an estimated semi-major axis of 3.4 AU and an eccentricity of 0.6 (Hatzes et 
al. 2000), orbiting with a period of ~ 7 yr. This inner planet clearly can not be responsible 
for the formation of the structures in the e Eridani debris disk, as it would not have any 
gravitational effect outside the inner ~ 20 AU of the system (Moran et al. 2004), but it 
confirms the presence of a planetary system around the star. 

Direct detection of e Eridani planetary companions has so far proved unsuccessful. Mac- 
intosh et al. (2003, hereafter M03) carried out near-IR Keck deep Adaptive Optics (AO) 
imaging that excluded the presence of planets with mass larger than 5 Mj at spatial scales 
comparable to the radius of the dust structures in the debris disk (roughly in a 40" search 
radius). Profntt et al. (2004, hereafter P04), while searching for disk scattered light in the 
optical with the Hubble Space Telescope (HST), compiled a list of ~ 60 objects between 
12.5" and 58" in a region E from e Eridani, none of which, in absence of a second epoch 
observation to measure proper motion, could be identified as a candidate substellar mass 
companion. 

The InfraRed Array Camera (IRAC, Fazio et al. 2004) on-board the Spitzer Space Tele- 
scope (Werner et al. 2004), thanks to its photometric system and unprecedented sensitivity 
in 4 photometric bands at 3.6, 4.5, 5.8 and 8.0 jum, can uniquely contribute to the search of 
substellar mass objects. Brown dwarfs and Jupiter-size planets are characterized by strong 
molecular absorptions features (CH 4 , NH 3 , H 2 0) in the mid-IR wavelength range. The IRAC 
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passband centered at 3.6 /xm, in particular, is very sensitive to a deep methane absorption 
feature at 3.3 jum, while the band at 4.5 /zm is largely free of molecular opacities. As a 
consequence, T-dwarfs and gas giants have nearly unique colors in the IRAC bands, and can 
be separated from other classes of astronomical sources in color-color and color-magnitudes 
diagrams, based on their IRAC photometry alone (Patten et al. 2004, 2006). Based on mod- 
els by Burrows et al. (2003), a 5 Mj planet at the distance of e Eridani will have a flux 
density of 1.3 mJy, and will be easily detected by IRAC if the artifacts due to the presence 
of the bright star can be mitigated. 

To take advantage of these unique characteristics of IRAC, we carried out deep imaging 
of e Eridani in all IRAC bands, as part of the "Fabulous Four debris disk star" Spitzer 
Guaranteed Time Observation (GTO) program (Program ID 90). This program consists of 
the observation of Vega, Fomalhaut, (5 Pictoris and e Eridani with all Spitzer instruments, in 
order to investigate the structure of their debris disks and to search for substellar companions. 
Results concerning the structure of the disks of Fomalhaut and Vega as seen by the Multiband 
Imaging Photometer for Spitzer (MIPS, Rieke et al. 2004) and the InfraRed Spectrograph 
(IRS, Houck et al. 2004) have already been published (Stapelfeldt et al. 2004; Su et al. 2005). 
A paper based on IRS spectroscopy and MIPS imaging of e Eridani disk is in preparation 
(Backman et al. 2006). We present here the result of the search for substellar companions 
around e Eridani made with IRAC. In Section 2 we describe the IRAC observations and the 
data reduction, with particular emphasis on the technique adopted to remove the e Eridani 
starlight diffracted on the focal plane, to allow detection of faint sources as close as possible 
to the star. In Section 3 we discuss our detection limits for disk extended emission and for 
low-mass companions. In Section 4 we describe the results of our companion search and 
compare them with previous searches. Finally, we present our conclusions in Section 5. 



2. Data acquisition and reduction 

Due to its proximity to the Sun, e Eridani is a very bright star in the mid-IR, with a 
2MASS K magnitude of 1.76 (Cutri et al. 2003). Faint companions and diffuse circumstellar 
emission can be detected only after a careful subtraction of the bright central point source. 
The Spitzer Space telescope is specially suited for this task, thanks to the exceptional stability 
of its optics and pointing system, that allows a precise measurement and reproducibility of its 
instrumental Point Spread Function (PSF). The IRAC observations for the whole "Fabulous 
Four" program were designed to maximize the ability of subtracting the stellar PSF. 

Each star was observed in two epochs, typically one month apart, to take advantage of 
the different roll angle of the spacecraft, which is determined by the relative position of the 
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spacecraft with the Sun. This strategy allows to position the main features of the PSF (the 
diffraction spikes created by the tripod supporting the secondary mirror and a few artifacts 
created by the electronics), which are fixed with respect to the detector pixel grid, in different 
positions on the sky. This allows a complete coverage of the sky around each of the stars. 

e Eridani was observed on January 9, 2004 (AOR 4876032) and on February 17, 2004 
(AOR 4876288). The roll angle offset between the two epochs is of 20.15°, clockwise. Each 
observation consisted in a sequence of 12 s Full Array frames (10.4 s integration time), 
dithered on the IRAC arrays on a 36 Position Reuleaux triangle, using the small dither 
scale, with each position repeated 9 times. The total integration time was 3379.6 s (over 56 
minutes) on-source, for each IRAC band. The source was set at the center of the array, where 
the PSF is cleaner from artifacts. The other "Fabulous Four" stars were observed with the 
same technique, and with a total integration time scaled according to the brightness of the 
source. A standard star (e Indi) was also observed, with the same total integration time as 
e Eridani, to contribute to the measurement of the instrumental PSF. The total integration 
time was estimated on the basis of the sensitivity required to detect a Jupiter mass planet 
in the vicinity of e Eridani. The individual frame exposure time (10.4 s) was chosen because 
it allows to get a long total integration time divided in a manageable number of individual 
frames, while still avoiding excessive saturation of the star on the detector. Shorter frame 
times are available, but could not be used because it would have been very inefficient to 
reach the required total integration time and the large field of view necessary to search for 
widely separated low mass companions. The total field of view imaged in each band was 
5.78', slightly larger than the IRAC field of view of 5.21', due to the dithering pattern. 

Basic data reduction and calibration was done with the Spitzer Science Center (SSC) 
pipeline, version S10.0.1 (first epoch) and S10.5.0 (second epoch). The IRAC pipeline returns 
the individual exposures calibrated in physical units of surface brightness. This calibration is 
valid only for point sources, as it was obtained by matching the flux of standard stars modeled 
by Cohen et al. (2003) within an aperture of 12.2". Note that point source photometry with 
different sets of apertures and sky annuli and the photometry of extended sources requires the 
use of appropriate aperture corrections. These aperture corrections and the other parameters 
used in the absolute photometric calibration (including the "FLUXCONV" factors used to 
convert the raw data in DN/s into the science calibrated data in MJy/sr) are listed in Table 1. 

For each epoch, separately, we have created a mosaic combining the individual frames, 
using the SSC Mosaicer MOPEX on a final grid with a pixel scale of 0.4"/pix (a factor 3 
smaller than the original IRAC pixel scale, of ^1.22"/pix) to leverage the high coverage for 
a better sampling of the PSF. We have chosen to maintain the same orientation in the final 
mosaic as the individual frames, to have the PSF features oriented with the same angle in 
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the mosaics of the two epochs. Cosmic rays and other outliers have been removed using the 
MOPEX temporal outlier module. This procedure was repeated identically for all "Fabulous 
Four" sources, including the standard star, in all epochs. 

2.1. PSF subtraction 

Figure 1 shows the "Fabulous Four" program PSF we have derived by combining the 
mosaics of Vega (2 epochs) and e Indi (1 epoch). Note that Vega itself is a star surrounded 
by a debris disk, but the infrared excess in IRAC bands is too small to be detectable by 
IRAC (Su et al. 2005). Vega can thus be considered as a point source for the purpose 
of building our PSF. The three mosaics have been shifted by matching the position of the 
unsaturated diffraction spikes, rescaled according to the relative brightness of Vega and e Indi 
and then coadded. The final PSF is clean of background stars, that have been a-clipped in 
the coadding phase. 

A special version of the PSFs we have generated for this project, that combines all the 
stars of the "Fabulous Four" program with the exclusion of (3 Pictoris (whose debris disk is 
indeed extended at IRAC wavelengths), is available at the SSC web site 1 for different pixel 
scales. These public PSFs also include an unsaturated core derived from the observation of 
faint standard stars as part of the IRAC photometric calibration project. A paper describing 
the construction of these PSFs and their main characteristics is in preparation (Marengo et 
al. 2006). 

It is clear from Figure 1 that the IRAC PSF has a lot of structures and artifacts that 
are the main limiting factor for observing faint emission around bright sources. These arti- 
facts cover the whole area of the IRAC arrays, and can be divided in two categories. The 6 
diffraction spikes (each of them composed by two initially braided and then diverging indi- 
vidual spikes), the smooth extended PSF tails, the diffraction rings and the PSF "ghosts" 
("filter ghosts", small "ring-like" and "cross-like" structures at the left of the PSF peak 
at 3.6 and 5.8 /im and the right at 4.5 and 8.0 /im, and "beam splitter ghosts", fainter 
structures below the 5.8 and 8.0 /im PSF peaks) are part of the "optical PSF". These struc- 
tures are linear in intensity with the fluency of the source (product of its surface brightness 
with the integration time) and can thus be easily subtracted even with a PSF derived from 
stars of different brightness than the source, as long as the PSF stars are observed with 
the same dither pattern and in the same position on the IRAC detectors (the position of 
the ghosts and the orientation of the diffraction spikes does depend on the array position). 
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-7- 



The negative columns ("column pull-down", white area in Figure 1) and the positive fencing 
pattern ("muxbleed", visible in Figure 9 after PSF subtraction) in the 3.6 and 4.5 /iin bands, 
the large positive "crosses" (vertical and horizontal "banding") and the saturated multiple 
"echos" on the right of the saturated core ("bandwidth effect") in the 5.8 and 8.0 jum bands 
are instead "electronic artifacts" generated by the detector and the readout electronics and 
artifacts resulting from internal scattering of IR light within the detectors. The column pull- 
down also affects the 5.8 and 8.0 /im bands, but can be hidden below the stronger banding. 
The electronic artifacts are not linear with the source fluency and can only be subtracted 
using a PSF obtained from stars with the same signal as the source. 

The PSF that we have derived from Vega and e Indi is only effective for subtracting 
the optical structures of the PSF. To eliminate the electronic artifacts it is necessary to 
use an alternative approach, based on the observation of e Eridani in two epochs, with two 
different roll angles. By using the second epoch as a PSF model for the first epoch (a PSF 
with exactly the same fluency of the source), we can effectively subtract the optical PSF 
and the electronic artifacts, while preserving the background sources that are rotated in 
the two epochs with respect to the PSF. In a two epoch subtracted frame the background 
sources from one epoch will be positive, and the ones from the other epoch will be negative. 
The trade-off for this technique is that for a very crowded field there is the risk of aliasing 
"positive" with "negative" sources when they superpose by chance, due to the rotation of 
the field. 

To have the best possible coverage of the sources in the e Eridani field, we have applied 
both PSF subtraction techniques in all IRAC bands. For each band we have thus produced 
4 PSF-subtracted frames: (1) first epoch subtracted with the PSF; (2) second epoch sub- 
tracted with the PSF; (3) first epoch subtracted with the second epoch and (4) second epoch 
subtracted with the first epoch (3 and 4 are of course the same, except for the orientation 
on the sky and the sign). 

Figure 2 shows an example of the two PSF subtraction techniques (right panels) com- 
pared with the original images before PSF subtraction (left panels). The top right panel 
shows the epoch 2 image subtracted with the PSF: note the remaining "cross" due to im- 
perfect subtraction of the "electronic" artifacts. The bottom right panel shows the second 
epoch image subtracted with the first: note the better subtraction of the "electronic" ar- 
tifacts, but the occurrence of aliasing for some of the field sources. Both PSF subtraction 
methods are able to remove most of the e Eridani light scattered by the PSF, and bring forth 
the numerous fainter sources in the field. Limits to this technique due to residual noise after 
the PSF subtraction are described in the following section. 

The PSF subtraction technique described here allows to measure the brightness of e Eri- 
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dani with considerable precision, by determining the scaling factor between the e Eridani 
image and the PSF (which is normalized, by construction, as the calibrated image of Vega). 
This scaling factor is determined by the fitting routine better to than 1%, which translates 
into an accuracy of ~0.01 mags. The measured magnitudes, in each epoch, are listed in 
Table 2 (the conversion from magnitudes to fluxes uses the IRAC Vega fluxes also reported 
in the IRAC Data Handbook 2004). Note that there is no photometric variation between 
the two epochs, with the exception of ~ 2.7% change in the 5.8 fj,m flux (which is however 
within 3a of the estimated photometric error). 



3. Sensitivity limits for extended and point sources 

The PSF subtraction methods described in the previous section are able to remove 
most of the e Eridani light scattered on the IRAC arrays by the Spitzer optics and the 
detector electronics. The subtraction, however, leaves a number of artifacts and residual 
noise that limits our capability to search for faint companions and extended structures. 
These residuals are especially strong in close proximity to the star, precluding the search for 
planetary companions in the inner e Eridani system. 

Figure 3 shows the radial profiles of e Eridani before (thin line) and after (thick line) 
PSF subtraction, along a wedge with Position Angle 120° and width 40° (positioned to avoid 
diffraction spikes and electronic artifacts). Note the presence of a background pedestal, in 
part due to an imperfect subtraction of the dark frames, and in part due to an electronic 
pedestal produced when a bright source is incident on the detector. This background pedestal 
is particularly strong in IRAC images containing bright sources, as it is the case with e Eri- 
dani, and is larger at longer wavelengths. The pedestal level in our observations was ~ 0.1, 
0.2, 2.5 and 6.0 MJy/sr at 3.6, 4.5, 5.8 and 8.0 fim respectively, and was successfully removed 
from the PSF-subtracted images (but not from the profiles shown in Figure 3) during PSF 
subtraction by adding a constant value to the scaled PSF, even though a small residual slope 
was left in the final 5.8 //m images. 

The IRAC arrays become severely non linear when the signal exceeds the ~ 30, 000 - 
50,000 DN range (see Table 1). Above this limits the image is rapidly saturated. This effect 
can be observed in the radial profiles as a surface brightness depression in the inner 5-7" (16- 
22 AU), where all signal is lost, making the core of the star looking like a "donut". Outside 
this radius the PSF profiles rapidly fall, even though the residuals after PSF subtractions 
are still high within ~ 14" from the star (~ 45 AU). These residuals are mainly due to 
the high photon noise of the bright source, the buildup of strong latents during the long 
total integration, and to the intrinsic limitations in the PSF construction method: even 
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small changes in the dither pattern, combined with the undersampling of the IRAC PSF 
(especially severe at 3.6 and 4.5 microns), can cause small but significant deviations in the 
PSF profiles, that result in elevated residual noise where the PSF signal is stronger. Outside 
a ~ 14" radius the PSF subtraction is however very effective to lower the PSF signal to the 
photon noise level. 

3.1. Limits on the debris disk detection 

P04 have used the STIS CCD camera on the Hubble Space Telescope to search for an 
optical counterpart of the sub millimeter disk in scattered light. The result of their search was 
inconclusive, but they set an upper limit of 25 mag/ arcsec 2 for the optical surface brightness 
of the dust at the distance of 55 AU from the star (17", roughly the distance where the 
millimetric disk emission peaks), within the "clear aperture" wide spectral bandpass of the 
STIS camera (from 0.2 to 1.02 /im). At the nominal wavelength of 0.7 /im, this limit 
corresponds to a surface brightness of 0.011 MJy/sr. 

Given that the spectral energy distribution of a K star like e Eridani peaks in the 
near-IR, the intensity of the light scattered from the disk rapidly decreases at wavelengths 
longer than 1 /im, at least until the thermal emission of small grains starts to be significant 
(for A > 10 yum). This means that at IRAC wavelengths the total disk surface brightness 
will be very small. P04 have developed a model of the disk surface brightness based on a 
previous model by Li et al. (2003), fitting the 850 /im and IRAS emission and including the 
contribution from scattered optical light. This model, that assumes the presence of highly 
porous particles, and a rather flat grain size distribution from 1 fim to 1 cm, predicts that 
the surface brightness of the disk at 3.6 //m and at the distance of 55 AU would be less than 
0.002 MJy/sr, and even lower in the other IRAC bands (see Figure 11 in P04). This surface 
brightness is too small to be detected in our PSF-subtracted images. 

Figure 4 shows the sensitivity of our images for extended sources, after subtraction of 
the e Eridani starlight. We initially computed the RMS noise in a running boxcar of 4" 
size for the PSF and two-epochs subtracted images, obtaining 4 maps for each IRAC band. 
The 4 maps differ in that the residuals of the PSF subtraction, and thus the residual noise, 
depend on the PSF subtraction method and on the rotation of the field with respect to the 
array coordinates. Since we are interested in the best case sensitivity achievable combining 
both epochs and both PSF subtraction methods, we have constructed a final map by taking 
the minimum RMS value, in each pixel, of the individual maps. We have then corrected 
the minimum RMS maps for the "infinite aperture correction" necessary to recalibrate the 
maps for the case of extended sources (listed in Table 1). 
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The contours show that the residual RMS noise increases closer to the star, as expected, 
and along the main PSF artifacts (diffraction spikes and electronic artifacts). The large oval 
shape at the left from the star at 4.5 /mi is enhanced noise due to a pupil ghost. At the 
distance of 55 AU from the star the sensitivity is greater than 1 MJy/sr in all channels, 
making impossible to test the accuracy of the P04 model. Even as far as 1' from the star, 
the maximum sensitivity (0.01 MJy/sr at 3.6 and 4.5 /mi and 0.05 MJy/sr at 5.8 and 8.0 /mi) 
is still not sufficient to detect the scattered light from the disk. 

The P04 model predicts a total flux of 0.12, 0.14, 0.19 and 0.40 mJy at 3.6, 4.5, 5.8 
and 8.0 /mi. Sheret et al. (2004), using a simple model of a thin ring with solid silicate 
and organic grains with a size distribution similar to the interstellar medium (Mathis et al. 
1977), but ranging from 1.75 /mi to 5 m, found a larger flux of 1.7, 1.5, 0.97 and 0.53 mJy 
in the IRAC bands. In both cases the predicted fluxes are much smaller than the accuracy 
of our IRAC photometry (see Table 2), and the expected infrared excess from the disk over 
the photospheric emission cannot be detected by our measurement. 

3.2. Limits on low mass companions detection 

Models shows that, at least in the case of a faint debris disk like the one around e Eridani, 
the detection of scattered light from the circumstellar dust is beyond the sensitivity limits 
of our PSF-subtracted images. The main goal of IRAC imaging, however, is the search for 
faint companions with sub-stellar mass. From the profiles shown in Figure 3 it is clear that 
we have no point source sensitivity closer than ~ 14" (45 AU) from the star, because of 
saturation, nonlinearity and high photon noise preventing usable PSF subtraction very close 
to the star. This rules out the possibility to search for the radial velocity planet at 3.2 AU, 
which would be well within the saturated area in our images. A different story, however, 
concerns the possibility to detect the planet whose presence at a distance between 40 and 
60 AU (12" - 20") is predicted by some dynamic models of the observed clumps in the 
debris disk. Even though such planet has been unsuccessfully searched by M03, the superior 
sensitivity of IRAC and the photometric system more suited for the observation of low mass 
objects may give us the chance of detecting this body, if it exists and if its mass is above 
our detection limits. The wide field of view of the IRAC maps (covering an area with size 
over 1,200 AU around the star) also allows us to search for widely separated planets that 
may have been ejected from the system, due to planet-planet interactions (Ford et al. 2001) 
in case the total protoplanetary disk mass was not sufficient to dampen the eccentricity of 
planets forming in the system. 

Since the discovery of Brown Dwarfs (Becklin & Zuckerman 1988; Nakajima et al. 2005), 



- 11 - 



several groups have started to compute detailed model atmospheres of substellar mass ob- 
jects to predict what would be their emission spectra and luminosity. More recently, these 
models became very sophisticated, to include large lists of molecular and alkali opacities, 
dust formation and settling, improved treatment of line pressure broadening and the effects 
of water clouds and ice formation, enabling them to compute reliable spectra extending in 
the IR domain. Figure 5 shows a family of such models, computed by Burrows et al. (2003) 
for 1 Gyr old substellar objects with mass from 1 to 25 Mj. The model spectra are shown 
in the wavelength range covered by the IRAC photometric system, and have been rescaled 
for the distance of e Eridani. As mentioned in section 1, the spectrum of substellar objects 
in the IRAC bands is modulated by deep CH 4 and H 2 absorption features. Methane ab- 
sorption is particularly strong in the 3.6 /im IRAC band, where it can depress the flux by 
a few orders of magnitude below the continuum, with increasing strength with decreasing 
temperature (and thus decreasing mass, for a fixed age). A second broad methane feature 
is present within the 8.0 /im band, while water absorption affects the flux around 5.8 [im. 
The 4.5 ^m IRAC band, however, is relatively free of molecular absorption, and thus offers 
the best chances of detection for molecular rich, cool substellar objects. 

We have estimated the magnitude and colors of the atmosphere models presented in 
Burrows et al. (2003) in the IRAC photometric system by integrating each model with the 
IRAC bandpasses, normalized by the flux of Vega: 



m (band) = -2.5 log 



Land Fl Ve9a) (A) Rv(\) A d\ 



(1) 



where F U (X) is the model spectrum, Fl Vega \X) is the model spectrum of Vega, derived by 
Robert Kurucz for Martin Cohen (this is the same model used for the absolute photometric 
calibration of IRAC), and -R^(A) is the IRAC Spectral Response Curve from the IRAC Spitzer 
Science Center web site 2 (August 9, 2004 version). 

The predicted magnitudes are shown in Figure 6, and compared with the 5a point 
source sensitivity of our PSF-subtracted images along the same radial cut used in Figure 3. 
We have derived the point source sensitivity from the maps in Figure 4 by converting the 
RMS noise into Noise Equivalent Flux Density (NEFD) inside an aperture of 2 IRAC pixels 
radius (2.4"): 



2 http://ssc. spitzer. caltech.edu/irac/spectral_response. html 
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rn 



^(x,y) = -2.5 log 



5-NEFD(x,y) 

F [Vega) 



(2) 



and 



NEFD(x, y) = 2.3504 • 10" 5 • r A ■ ^(pixscale) 2 ■ acorr{rA) ■ RMS(x, y) 



(3) 



where is the point source sensitivity in magnitudes at the (x, y) coordinates in the 

PSF-subtracted images, RMS(x, y) is the sensitivity map shown in Figure 4, F^yega) is the 
flux of Vega in the IRAC bandpasses reported in the IRAC Data Handbook (2004), and 
the other quantities are geometrical factors necessary to convert surface brightness into flux, 
and include the aperture radius r^, the pixel scale of the PSF-subtracted images (0.4"/pix) 
and the aperture correction acorrijA) for 2 IRAC pixel apertures, also from the IRAC Data 
Handbook (2004), listed in Table 1. 

Figure 6 shows that indeed the IRAC band most sensitive to detect cool companions is 
the one at 4.5 jj,m. Outside the area where the image is contaminated by saturation, non- 
linearities and photon noise from the central bright star (cross-hatched in the figure), our 
PSF-subtracted images are sensitive to Jupiter mass planets. Just inside the inner edge of 
the ring we are sensitive to planets with a mass of ~ 1 Mj, and below 1 Mj in the area outside 
the debris disk. The strong methane absorption in the 3.6 /xm band increases considerably 
the minimum detectable mass at this wavelength, so much that in the region inside the 
disk where dynamic models suggest the presence of planets gravitationally perturbing the 
disk, we are only sensitive to bodies with mass larger than 7 Mj. The sensitivity in the 
two remaining IRAC bands (5.8 and 8.0 /mi) is intermediate between these two extreme 
cases. Note that, with the exception of the 3.6 /im band, outside a radius of 45 AU our 
PSF-subtracted images are at least a factor 2 more sensitive in terms of mass than the Keck 
K band images in M03, and at least 5 times more sensitive outside the submillimeter ring. 



Even though e Eridani is relatively far from the galactic plane, the field around the star 
shows a large number of background sources. The challenge is determining if any of them 
can be a gravitationally bound low mass companion, instead of a faint background star or 
galaxy, or an artifact from the PSF subtraction. The similar search made by P04 in the 
optical concluded that most of the detection were in fact background galaxies. Given the 
sensitivity of IRAC to red extragalactic objects, we should expect a similar result. 



4. Companion search 
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As mentioned above, the colors of substellar objects in the IRAC photometric system 
are very peculiar. The colors of brown dwarfs and giant planets at the distance of e Eridani 
can be estimated from Burrows et al. (2003) models, as described in section 3.2. For all 
background sources having good photometry in at least two or three IRAC bands, it is 
thus possible to investigate their nature on the basis of color-color and color-magnitude 
diagrams. Planets below 5 Mj, however, may be detected in only one IRAC channel, at 
4.5 /iin, because of the lower sensitivity at 5.8 and 8.0 ^m, and the deep methane absorption 
feature depressing their 3.6 ^m flux below our sensitivity level. Given that our two epochs 
are too close in time to allow for proper motion measurements, for these "3.6 ^m dropout" 
sources a new observation at a later epoch is needed to check if they are comoving with the 
e Eridani system. 



4.1. Source detection and photometry 

The PSF subtraction procedure leaves behind a large number of point-source-like arti- 
facts that can easily confuse the detection algorithms available in most photometry packages. 
In particular, analysis based on the shape of the point sources by Gaussian or PSF fitting of 
all structures above a certain threshold would reject a large number of valid sources whose 
shape has been distorted by the proximity with the PSF artifacts, or by a nearby "aliasing" 
source in the two epoch subtracted frames. Given our requirement to be inclusive, not to 
miss the proverbial "needle in the hay stack" that would be a planet orbiting e Eridani, we 
have relaxed our requirements starting with a list including any detection with S/N > 5 
within an aperture of 2.4" radius. We have then selected the real from the spurious sources 
by visually inspecting each detection, and by comparing the two epochs and the results from 
the two different PSF subtraction methods described in section 2.1. With this procedure we 
have obtained a list of 467 sources with magnitudes between 12 and 20, reliably detected in 
at least one epoch at 4.5 ^m within the IRAC field of view. 

Fazio et al. (2004b) shows that, outside the plane of the Galaxy, most of the sources 
detected by IRAC with a magnitude > 13 at 4.5 fim are background galaxies. Table 1 in 
that paper shows that the total number of background galaxies with 4.5 /im magnitudes 
between 15 and 20 within the area of our maps is ~ 920. This is consistent with our 
number of detections, taking into account the large area at the center of our PSF-subtracted 
images where our sensitivity limits is well below the 20th magnitude because of the high 
PSF subtraction residuals. 

For source detection and photometry we have used the package "PHOTVIS" developed 
by Rob Gutermuth (Gutermuth et al. 2004), based on the "APER" aperture photometry 
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package part of the IDL "ASTROLIB" library. The size of the aperture was set to 2.4", 
with a sky annulus with inner and outer radii of 2.4" and 7.2". To convert the mosaic data 
(in MJy/sr surface brightness units) into magnitudes we have derived zero point magnitudes 
based on the Vega fluxes adopted for the camera's absolute calibration, listed in Table 1. 

The uncertainty of each photometric measurement was estimated as the sum in quadra- 
ture of the scatter in the sky values, the uncertainty in the mean sky brightness and the 
random photon noise on source (estimated from the source counts converted in electrons). 
Note that this error analysis tends to underestimate systematic errors due to imperfect sub- 
tractions of the PSF and PSF artifacts that alter the sky level and introduce non-Gaussian 
noise. 

We have then merged the individual source lists in each band and PSF subtraction 
method. In each IRAC band we have averaged the individual photometry in the two epochs 
and PSF subtractions, weighted on their S/N ratio. Typical uncertainties of the final pho- 
tometry are ~ 0.05 mag at 3.6 and 4.5 /im, and ~ 0.1 mag at 5.8 and 8.0 /im. Sources close to 
e Eridani or along PSF subtraction artifacts may have a larger uncertainty in one of the two 
epochs, and depending on the PSF subtraction method. We have checked the consistency 
of the photometry of sources detected in both epochs, and found that it is generally within 
the error margin, except when a PSF artifact or source aliasing is present. In these case, we 
assume an uncertainty of ~ 0.2 mag, representative of the maximum spread observed in the 
photometry of sources detected in both epochs in at least two IRAC bands. 



4.2. Color-color and color-magnitude diagrams 

The top panel of Figure 7 shows the [3.6]-[4.5] vs. [4.5]-[5.8] color-color diagram of the 
sources detected in all three bands. Note that many of the sources are clumped around 
zero colors, indicating that they are regular background stars. A significant number of 
sources, however, have one or both colors red: this is expected for extragalactic sources 
having significant PAH emission in the longer wavelength IRAC bands, or red-shifted galaxies 
(Huang et al. 2004; Barmby et al. 2004). The bottom panel shows the [3.6]-[4.5] color vs. the 
4.5 /im absolute magnitude that the sources detected in both 3.6 and 4.5 /im bands would 
have if they were at the same distance (3.22 pc) of e Eridani. 

The contours plotted on the figure enclose the locations where sources with the colors of 
field M, L and T dwarfs are located in the diagram. These contours have been derived from 
data collected as part of the IRAC Guaranteed Time observations of M, L and T dwarfs 
(Patten et al. 2004, 2006), using a variation of the weighted ^-Nearest Neighbor Method 
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(kNN, Fix & Hodges 1951) adapted for astronomical applications. The method, commonly 
used in diagnostic medical studies is a multi-variate analysis very effective in the automated 
classification of items in different groups. In this work we have used k = 3 (3rd nearest 
neighbor) and a Gaussian kernel to bias the metric towards the 1st nearest neighbor. Sources 
inside each region differ from the class templates (the M, L, and T field stars observed by 
Patten et al. 2006) by less than la error in photometry or color. A detailed description of 
this method will be given in a separate paper (Marengo & Sanchez 2006, in preparation). 

The top panel of Figure 7 shows that many of the sources detected around e Eridani 
have colors compatible with the colors observed for field M, L and T dwarfs (Patten et al. 
2006). The bottom panel of the same figure, however, shows that all the sources detected 
at least at 3.6 and 4.5 /im are much fainter than all the field M, L and T dwarfs in the 
Patten et al. (2006) sample. Although the observed field M, L and T dwarfs have a range 
of masses and ages which may not be appropriate for e Eridani, lower mass brown dwarfs 
and planets at a given [3.6]-[4.5] color cannot be much fainter than the M 4 5 range shown 
for the field objects. The reason is that the radii of ~ 1 Gyr brown dwarfs and planets span 
a relatively small range of values, and are relatively constant independent of age and mass 
(Burrows et al. 1997; Chabrier & Baraffe 1997). Hence their luminosity depends primarily 
on temperature, and consequently color. 

This is demonstrated in Figure 8 by the Burrows et al. (2003) models of 1 Gyr old brown 
dwarfs and giant planets. The large area covered by the /cNN contours is determined by the 
large systematic error we have assumed for the models, approximately as large as the color 
separation between the models of individual mass. This is a prudent assumption, in our 
requirement to be inclusive, given that these models have never been directly tested towards 
real extrasolar planets. Even so, none of the detected sources possess the right combination 
of colors and magnitudes to be a candidate planet or T dwarf, as all the planet models have 
a [3.6]-[4.5] color much redder, and a [4.5]-[5.8] color much bluer (due to H 2 absorption at 
5.8 fim) than the sources detected in at least two IRAC bands. 

4.3. The 3.6 /im dropout sources 

Figure 8 shows that a planet with 5 Mj would have a 4.5 /xm absolute magnitude of 
~ 15.3 and a [3.6]- [4.5] color of 4.9. This means that the 3.6 ^m magnitude of that planet 
at the distance of e Eridani would be 20.2, which is below our detection limit. Such a 
planet can be detected in the 4.5 /im images, but would drop out from our 3.6 /im images. 
Bodies with mass below < 5 Mj cannot be found on the basis of their color with the A;NN 
method described in the previous section, and their identification should rely on proper 
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motion measurements. 

The time lag between our two epochs, however, is too short to yield a measurable shift 
of e Eridani and any companion. The proper motion of e Eridani is of -0.977"/yr in R.A. and 
0.018"/yr in Dec (Perryman et al. 1997), which means that during the 39 days between the 
two epochs the star moved by ~ 0.1". This is a factor 2 less than our best case centroiding 
accuracy (~ 0.2" at 4.5 fim, where the IRAC pixel scale has the optimal PSF sampling). 
Comparison with other catalogs is not possible, with the exception of the few sources in 
common with the P04 and M03 surveys, as described in the next section. 

Table 3 lists all the sources, within a radius of 600 AU (approximately the area covered 
by our images) that have a 5cr detection at 4.5 jum, but that are not detected at the 5<r level 
at 3.6 fxm. Most of these sources do have, however, a marginal 3.6 yum detection, and their 
colors (even allowing for the large uncertainties in their 3.6 yum photometry) suggest them to 
be red background galaxies. Only three sources are completely undetected at 3.6 ^m. Their 
4.5 //m fluxes and [3.6]-[4.5] colors (estimated from the 3.6 /xm sensitivity limit in Figure 4) 
are compatible with the colors and magnitudes of planets with masses < 3 Mj (sources 5 
and 28 in Table 3) and < 2 Mj (source 6 in Table 3), but not with faint red-shifted galaxies 
(which typically have [3.6]-[4.5] < 0.5). Huang et al. (personal communication) measured 
that out of 7.3 ■ 10 4 sources detected by IRAC in an area of 2° x 10' with a 3 hours exposure 
only 6 have a 4.5 /mi magnitude between 17.5 and 20, and [3.6]- [4.5] > 1.5. It is thus very 
unlikely to have 3 such sources within a single IRAC field of view. A careful inspection of 
the individual 4.5 y«m PSF-subtracted images however suggests these sources to be multiple 
(and possibly extended, at least in the case of source 5), which may indicate that their colors 
may have been altered by chance superposition with higher than average PSF subtraction 
artifacts. The real nature of these three detections cannot be determined with the current 
data, and a new observation to confirm that they are real, and to measure their proper 
motion, is necessary to determine if they are physically associated to e Eridani. 



4.4. Comparison with Keck AO and HST searches 

Even though none of our detected sources (with the possible exception of the three 
3.6 ^m "dropout sources" described in the previous section), do not have the colors or 
magnitudes to be low mass companions of e Eridani, it is however useful to compare the 
result of our survey with the previous M03 search. Figure 10 shows the location, in our 
4.5 yum images, of the 10 objects identified by M03 in their K band Keck AO maps within an 
area roughly as 40" (~ 130 AU) from e Eridani. We detected four of these objects (numbered 
3, 5, 6 and 10) in at least 2 IRAC images. We could not detect sources 1, 2, 4, 7, 8 and 



-17- 



9, which are below the PSF-subtraction residual noise in their position. We have however 
detected one source that was missed in the Keck search (inside the square box number 11). A 
careful examination of the Keck plates (Figure 2 in M03) indeed shows the presence of a faint 
source (below the Keck observations detection limit) in the location of source 11. Sources 
4, 6 and 11 have also been detected by P04: they have an extended sources appearance and 
are probably background galaxies. 

Table 4 summarizes the properties of the M03 sources. The "Keck Coordinates" columns 
indicate the position of each source at the time of the M03 first epoch (December 2001). Since 
Table 3 in M03 does not provide absolute coordinates for the sources, but only offsets from 
e Eridani, we had to derive their actual R.A and Dec by estimating the position of e Eridani in 
December 2001 from the ICRS 2000 (Perryman et al. 1997) coordinates and proper motion. 
Given that M03 had to keep e Eridani outside their field of view to limit PSF artifacts, 
the accuracy of the December 2001 coordinates is limited to 0.2", which is the precision 
with which M03 can establish the position of e Eridani relative to their images. We then 
estimated the coordinates of each source at the time of the first IRAC epoch (January 2004) 
by extrapolating the motion measured by M03 between December 2001 and their second 
epoch in August 2002. Note that for some of the sources, either because they were used 
as position reference, or because the second epoch was missing (source 3), the predicted 
2004 position cannot be derived. Given the high accuracy quoted by M03 in the relative 
shift measurements, the uncertainty of the predicted coordinates in the IRAC images is still 
dominated by the uncertainty in their 2001 absolute coordinate, of 0.2". 

Finally we determined the position of all sources detected within the Keck field of 
view by measuring their photo-center in the IRAC 4.5 /im images. The uncertainty in 
the measurement of the sources centroid is limited by the presence of PSF subtraction 
artifacts and second epoch source aliasing, but is consistent, between our two epochs, with 
an accuracy of 0.2". Within this limit, we didn't observe any shift in the source coordinates 
between the two IRAC epochs. The real uncertainty in the IRAC coordinates, however, is 
much larger than 0.2", due to systematic errors in the pointing of the spacecraft. These 
errors are normally corrected by a pointing refinement procedure, included in the Spitzer 
Space Center data reduction pipeline, that uses 2MASS stars (Cutri et al. 2003) to find an 
absolute reference frame for the IRAC images within the 0.15" 2MASS accuracy. In the 
case of e Eridani, however, this pointing refinement was not done because there isn't any 
2MASS source detected in the vicinity (because of the e Eridani glare saturating the 2MASS 
detectors). The headers of the IRAC BCD frames show a pointing accuracy error of ~ 1" at 
the time of our observations. 

By comparing the coordinates of all sources detected in both Keck and IRAC images, we 
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have measured the shift between December 2001 and January 2004. The values we obtained, 
listed in Table 4 are rather large, and appear to suffer a systematic shift, possibly due to the 
limitations in the Spitzer pointing accuracy described above. We have tried to correct this 
systematic shift by subtracting the average value of the the measured shift, of 0.13" in R.A. 
and 1.30" in Dec. We estimate the uncertainty of the final shifts to be of the order of ~ 1". 
A significant fraction of this uncertainty may be due to the different centroid that partially 
resolved extragalactic sources, if extended or multiple, may have in the K and 4.5 fim bands. 
Within this error, none of the sources show a measurable proper motion, and the corrected 
shifts listed in Table 4 are all well below the expected e Eridani proper motion between the 
Keck and IRAC observations, of more than 2" in R.A. We also did not find any significant 
shift between the coordinates of sources 4, 6 and 11 as measured in the IRAC 4.5 fim images, 
and in the STIS/HST P04 maps. 

This confirms that none of the sources detected in the common field of view with the 
Keck and P04 observations is a low mass companion of e Eridani. 

5. Conclusions 

As the number of exoplanets detected with radial velocity or eclipsing techniques, since 
the initial detection by Mayor & Queloz (1995), continue to grow, the payoff of directly 
detecting their thermal emissions cannot be underestimated. As of today, thermal radiation 
has been directly detected from only three planets, all of them eclipsing "hot Jupiters" 
whose light has been detected with Spitzer: TrES-1 with IRAC (Charbonneau et al. 2005) 
at 4.5 and 8.0 /an, HD209458b (Deming et al. 2005) with MIPS 24 ^m and HD189733b with 
IRS (Deming et al. 2006). These three detections allowed for the first time a comparison 
between theoretical models of hot Jupiters with real data. However, no direct detection of 
gas giants orbiting a nearby parent star at a distance comparable to our own Jupiter and 
Saturn has been achieved so far, as all searches conducted from the ground (M03) and from 
space (P04 and this one) have returned negative results. Models of giant planets orbiting 
systems analogous to the Solar System are thus still untested. 

While future space missions, such as the James Webb Space telescope and the Terrestrial 
Planet Finder are being designed with the goal of direct exoplanetary detection, this work 
shows that the Spitzer Space telescope with the IRAC camera already has the capability of 
detecting young gas giants with a mass as low as a few Mj if they are orbiting their parent 
star at a distance comparable to our own Solar System Kuiper belt. Planets ejected on very 
elongated orbits during the initial phases of planetary system formation can be detected 
with Spitzer/IRAC down to masses of 1 Mj. Our observations allow to set a limit of ~ 2 Mj 
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at the distance of ~ 60 AU (~ 18") from the star, and < lMj outside the debris disk ring 
(~ 35"~ 112 AU). 

Comparison with M03 confirms that within our search space we are able to detect all the 
sources found with the most sensitive ground based search outside a radius of ~ 130 AU and, 
in addition, one below the Keck sensitivity. Comparison with the position of the common 
sources between IRAC and Keck observations confirms that none of the detections match the 
proper motion required to be a physical companion of e Eridani, validating and extending 
M03 results. 

Based on the colors of the detected sources, we did not find any strong companion 
candidate. However, for a number of sources, listed in Table 3, we are missing a detection 
in the 3.6 fj,m IRAC band. Without a good [3.6]- [4.5] color we cannot determine the nature 
of these sources, and a second epoch is necessary to measure their proper motion. Given the 
positional accuracy in the PSF-subtracted IRAC frames of ~ 1", such a measurement will 
be feasible after a time lag of about 2 years from the current observations. 

We would like to thank Brian Patten for providing the IRAC photometry of the MLT 
dwarfs we used as templates for our /cNN method, Robert Gutermuth for providing the 
PHOTVIS software, Mayly C. Sanchez for helping to adapt the ANN method to this appli- 
cation, Jiasheng Huang for kindly providing unpublished data used to determine the statistics 
and IRAC colors of extragalactic sources and Michael Schuster for a careful reading of this 
manuscript. 

This work is based on observations made with the Spitzer Space Telescope, which is 
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Table 1. IRAC Photometric Calibration 



Item 


3.6 /iin 


IRAC band 
4.5 fj,m 5.8 fim 


8.0 fan 


Isophotal A [yum] 


3.458 


4.492 


5.661 


7.870 


FLUXCONV [(MJy/sr)/(DN/s)] 


0.1125 


0.1375 


0.5913 


0.2008 


Saturation limits [DN] 


30,000 


35,000 


45,000 


50,000 


Saturation limits [MJy/sr] 


320 


460 


2560 


970 


Infinite aperture corrections 


0.94 


0.94 


0.63 


0.69 


2.4" aperture corrections 11 


1.213 


1.234 


1.379 


1.584 


Zero point magnitudes 


19.46 


18.97 


18.38 


17.56 


F„(Vega) [Jy] 


277.5 


179.5 


116.5 


63.1 



a Based on IRAC Data Handbook (2004) 

b For sky annulus with 2.4" and 7.2" inner and outer radii 

Tor pixel size 0.4"/pix, include 2.4" aperture correction 
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Table 2. e Eridani IRAC magnitudes and fluxes 



January 9, 2004 February 17, 2004 
IRAC band A [/im] a mag b Flux [Jy] c mag b Flux [Jy] c 



[3.6] 


3.548 


1.59 


63.9±0.3 


1.59 


64.1±0.4 


[4.5] 


4.492 


1.66 


39.1±0.2 


1.65 


39.2±0.2 


[5.8] 


5.661 


1.61 


26.5±0.2 


1.64 


25.8±0.2 


[8.0] 


7.870 


1.60 


14.5±0.1 


1.59 


14.6±0.1 



a IRAC nominal wavelengths from the IRAC Data Handbook (2004) 
b Magnitudes accurate to ±0.01 mag 

c Using IRAC Vega fluxes from the IRAC Data Handbook (2004) 
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Table 3. Sources with marginal or no detection at 3.6 /im within 600 A.U. from e Eridani 



# 




RA 


Dec 


[3.6] 

L J 


[4.5] 

L J 


[3.6] -[4.5] 


dist (AU) 


1 


53 


.189347 


-9.444469 


18.0 ±0.6 


17.29 ±0.19 





.3 ±0.6 


512 


2 


53 


.194310 


-9.465338 


19.2 ±0.7 


18.65 ±0.07 





.5 ±0.7 


437 


3 


53 


.194881 


-9.465428 


19.0 ±0.7 


18.52 ±0.05 





.5 ±0.7 


430 


4 


53 


.201029 


-9.469290 


19.7 ±0.7 


19.59 ±0.15 





,1±0.7 


373 


5 a 


53. 


.205047 


-9.452265 


> 20.2 


17.30 ±0.17 




> 2.9 


312 


6 b 


53. 


.205817 


-9.486791 


> 20.5 


19.14 ±0.10 




> 1.4 


443 


7 


53. 


.206656 


-9.489741 


19.6 ±0.6 


19.32 ±0.11 





.3 ±0.6 


464 


8 


53. 


.208558 


-9.492060 


19.7 ±0.6 


18.96 ±0.06 





.7 ±0.6 


472 


9 


53. 


.213247 


-9.454328 


19.0 ±0.6 


18.66 ±0.07 





.3 ±0.6 


214 


10 


53. 


.213330 


-9.447145 


19.5 ±0.6 


18.43 ±0.07 


1 


,1±0.6 


245 


11 


53 


.213868 


-9.414921 


17.8 ±0.5 


17.70 ±0.08 





,1±0.5 


543 


12 


53. 


.219017 


-9.474964 


19.2 ±0.7 


18.42 ±0.07 





.8 ±0.7 


240 


13 


53 


.219744 


-9.481340 


19.4 ±0.8 


17.57 ±0.08 


1 


.8 ±0.8 


299 


14 


53 


.221378 


-9.448212 


18.7 ±0.8 


18.38 ±0.11 





.3 ± 0.8 


164 


15 


53. 


.223314 


-9.488778 


19.7 ±0.6 


18.50 ±0.07 


1 


.2 ±0.6 


366 


16 


53. 


.227599 


-9.478849 


19.0 ±0.6 


18.28 ±0.05 





.3 ±0.6 


242 


17 


53 


.229451 


-9.440029 


18.9 ±0.6 


17.85 ±0.05 


1 


,0±0.6 


213 


18 


53. 


.234574 


-9.413988 


18.6 ±0.4 


17.76 ±0.07 


0. 


.8 ±0.4 


516 


19 


53. 


.236470 


-9.497356 


20.7 ± 1.9 


19.13 ±0.09 


1 


.6 ± 1.9 


458 


20 


53 


.237040 


-9.504618 


19.3 ±0.7 


18.44 ±0.06 





.9 ±0.7 


542 


21 


53. 


.238661 


-9.440854 


19.4 ± 0.9 


18.46 ±0.06 





.9 ±0.9 


220 


22 


53. 


.240732 


-9.441868 


19.1 ±0.9 


18.63 ±0.09 





.5 ±0.9 


220 


23 


53 


.249281 


-9.498276 


19.1 ±0.9 


18.60 ±0.08 


0. 


.5 ±0.9 


509 


24 


53 


.253290 


-9.495511 


17.7 ±0.3 


17.54 ±0.03 





.2 ±0.3 


502 


25 


53. 


.258599 


-9.484475 


19.3 ±0.7 


18.48 ±0.05 





.8 ±0.7 


440 


26 


53 


.259360 


-9.471248 


19.1 ±0.6 


18. 97 ±0.09 





1±0.6 


359 


27 


53 


.268980 


-9.431036 


18.4 ±0.4 


17.70 ±0.03 


0. 


.7 ±0.4 


541 


28 c 


53 


.269965 


-9.487714 


> 20.2 


17.64 ±0.16 




> 2.6 


565 


29 


53. 


.272692 


-9.430506 


18.7 ±0.5 


18.34 ±0.07 





.4 ±0.5 


580 



a Multiplc or extended sources detected in all 4.5 /xm images. 

b Point source (possibly double) detected in second epoch 4.5 /im images. 

c Point source (possibly multiple) detected in second epoch 4.5 /im images. 



Table 4. Proper motion and photometry of M03 objects 



s, 

7r 


Keck Coordinates' 1 


IRAC Coordinates 


Predicted 13 


Measured b ' c 


Corrected 13 " 




Photometry 6 






R.A. 


Dec 


R.A. 


Dec 


Aa 


AS 


Aa 


AS 


Aa 


AS 


[Kl 


[3.6] 


[4.5] 


[5.8] 


[8 01 


e Eridani 


53.23216 


-9.45825 


53.23159 


-9.45824 


2.05 


-0.04 




















1 


53.22949 


-9.45431 


















17.3 








— 


2 


53.23341 


-9.45353 


















17.3 








— 


3 


53.23099 


-9.44600 


53.23085 


-9.44634 






0.50 


1.22 


0.37 


-0.08 


16.3 


15.4 


15.4 


15.5 


15.4 


4 f 


53.23960 


-9.45539 






-0.03 


0.01 










19.4 










5 


53.24291 


-9.45422 


53.24316 


-9.45455 


0.05 


0.28 


-0.90 


1.19 


-1.03 


-0.11 


20.7 


17.4 


16.5 


15.1 


15.7 


6 f 


53.24219 


-9.46369 


53.24212 


-9.46390 


0.09 


-0.03 


0.25 


0.76 


0.12 


-0.54 


20.2 


17.6 


17.1 






7 


53.22852 


-9.45369 






0.05 


0.16 










20.3 










8 


53.22343 


-9.45447 






0.08 


0.00 










20.1 










9 


53.23658 


-9.46461 






-0.10 


0.06 










20.8 










10 


53.22722 


-9.46772 


53.22703 


-9.46828 


0.06 


0.10 


0.68 


2.01 


0.55 


0.71 


19.3 


16.7 


16.2 




15.7 


ll f 






53.24301 


-9.45902 


















16.5 




15.3 



a Positions calculated for Keck first epoch, Dec 1, 2001 
b Shifts expressed in second of arc 
c Shifts statistical uncertainty is 0.2" 

d Shifts corrected by subtracting the average shift of 0.13" in R.A. and 1.30" in Dec 
e Average statistical error in IRAC photometry is 0.2 mag 
f Detected in P04 as extended structures 
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Fig. 1.— The IRAC PSF at 3.6 /zm (a), 4.5 /an (b), 5.8 /an (c) and 8.0 jum (d). The PSF 
is shown in inverted logarithmic color scale, from zero intensity (white) to saturation level 
(black). The size of each panel is the same as the IRAC field of view (~ 312"). The main 
artifacts are marked. The "beam splitter ghosts" at 5.8 and 8.0 /im are the faint "kinks" 
visible along the vertical banding where the arrows point. Note that there is column pull- 
down also at 5.8 and 8.0 /im, which is however not visible due to the predominant effect 
of the banding. Muxbleed is also present at 3.6 and 4.5 /tm, and can be seen after PSF 
subtraction in Figure 9. Note that the horizontal spikes at 5.8 and 8.0 jum are hidden below 
the banding. 
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Fig. 2. — e Eridani before (left) and after (right) PSF subtraction. In the color composite, 
blue is 3.6 /im, green is 4.5 /zm, yellow is 5.8 /im and red is 8.0 /im. Note the rotation of the 
field of view between epoch 1 (a) and epoch 2 (c). Panel (b) is epoch 2 subtracted with the 
PSF constructed from Vega and e Indi. Panel (d) is epoch 2 subtracted with epoch 1. The 
left panels are shown in logarithmic color scale, while the PSF-subtracted images are shown 
in linear scale from -0.5 MJy/sr to 1.0 MJy/sr. The total field of view is 334". 
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Fig. 3. — e Eridani (epoch 2) radial profiles before (thin line) and after (thick line) subtraction 
with the IRAC PSF in the four IRAC bands. The profile is a circular average within a wedge 
of 40° width, with a position angle of 120° from the vertical, clockwise (to avoid diffraction 
spikes and other artifacts). The background pedestal level has not been subtracted from the 
profiles presented in this figure (even though it is removed from the PSF subtracted images 
used in the rest of the analysis). 
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Fig. 4. — Sensitivity maps, in array coordinates, computed as the RMS noise of the PSF- 
subtracted images, estimated in a running boxcar of size 4". The contours show the RMS 
noise of 5, 1, 0.5, 0.1, 0.05 and 0.01 MJy/sr, with the highest noise level at the center. The 
dotted line indicates the direction of the cut used to plot the radial profiles in Figure 3 and 
6. 
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Fig. 5. — Model spectra of gas giants and Brown Dwarfs from Burrows et al. (2003), com- 
puted for 1, 2, 5, 7, 10, 15, 20 and 25 Mj 1 Gyr bodies. The predicted flux density has been 
rescaled for the distance of e Eridani. The main spectral features in the IRAC bands are 
produced by molecular absorption of CH4 (in the 3.6 /im and 8.0 /im bands) and H2O (in 
the 5.8 /im band). The band at 4.5 /mi is relatively free of absorption, and is where Jupiter 
size planets would appear brighter. 
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Fig. 6. — Detection limits for giant planets and Brown Dwarfs based on models computed 
by Burrows et al. (2003). The solid line shows the 5a point source sensitivity (within an 
aperture of 2.4" radius) derived from the RMS maps in Figure 4, along a cut with position 
angle 120° and 40° width. The cross-hatched region within a radius of 14"(45 A.U.) from the 
star should be excluded because of saturation and high residual noise after PSF subtraction. 
The dashed vertical lines show the region where dynamic models predict the presence of 
a resonant body (between 40 and 60 A.U.), while the dotted vertical line marks the outer 
radius of the sub-millimetric ring at 35" (112 A.U.). 
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Fig. 7. — Color-color and color- magnitude diagrams of all point sources detected in the field 
around e Eridani. The contours are the region selected with the A;NN method for M dwarfs 
(solid line), L dwarf (dotted line), early T dwarfs (T<5, dashed line) and late T dwarfs (T>5, 
dot-dashed line). Note that many sources have the right colors to be M, L or T dwarfs, but 
none of them has the correct brightness to be an M, L or T companion of e Eridani, according 
to the colors and magnitudes of field M, L and T dwarfs from Patten et al. (2006). 
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Fig. 8. — Color-color and color-magnitude diagrams of 850 Myr giant planet model tracks 
from Burrows et al. (2003), compared with the sources detected in the e Eridani field. Square 
points denote the colors and magnitudes of model planets with mass from 25 to 1 Mj. 
The solid line encloses the region in the diagrams where candidate planets may be found, 
according to our &NN method. None of the sources detected at both 3.6 and 4.5 fim has 
the correct colors and magnitudes predicted by Burrows et al. (2003) models. The arrows 
indicate the photometry and colors of the three 3.6 /iin "dropout sources" listed in Table 3. 
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Fig. 9. — PSF-subtracted IRAC 4.5 /im images of e Eridani background field: (a) epoch 1 
subtracted with the PSF, (b) epoch 2 subtracted with the PSF, (c) epoch 1 subtracted with 
epoch 2 and (d) epoch 2 subtracted with epoch 1. The figure uses an inverted linear color 
scale from background level to 0.1 MJy/sr. The inner circle marks the maximum radius of 
the sub-millimeter ring (112 AU, corresponding to 35" at the e Eridani distance). The outer 
circle marks the distance of 600 AU. Circle points mark the location of point sources for 
which only a marginal detection (< 5a) at 3.6 //m is available. The three sources with no 
detection at all at 3.6 jttm are indicated with square points. 



-36- 




Fig. 10. — IRAC 4.5 fim images of the area surrounding e Eridani that was searched by M03: 
(a) epoch 1 PSF subtracted; (b) epoch 2 PSF subtracted; (c) epoch 1 - epoch 2; (d) epoch 
2 - epoch 1. The large squares are the field of view of the individual Keck images (40" size). 
The circles marks the positions of the sources found in the M03 search extrapolated to their 
expected position at the epoch of the IRAC observations. The small box indicates a source 
found in our IRAC images that was missed in M03. The two little diamonds indicate the 
position of e Eridani at the time of the first Keck epoch (left diamond), and the IRAC (right 
diamond) observations. 



